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Abstract 

We adapt the inverse iteration method for symmetric matrices to some nonlinear 
PDE eigenvalue problems. In particular, for p € (l,oo) and a given domain 17 C 
M n , we analyze a scheme that allows us to approximate the smallest value the ratio 
\D , if}\ p dx/ Jq \ip\ p dx can assume for functions ijj that vanish on <917. The scheme in 
question also provides a natural way to approximate minimizing ip. Our analysis also 
extends in the limit as p —>• oo and thereby fashions a new approximation method for 
ground states of the infinity Laplacian. 


1 Introduction 

In this paper, we will use a generalization of the inverse iteration method for symmetric 
matrices to estimate solutions of some nonlinear PDE eigenvalue problems. The first problem 
we consider is as follows. For p e (1, oo) and a bounded domain 17 C M”, we define 

A ' : = u dS^ : ^ <p(n) ’^ 0 }' (i - i) 

Here W 0 1,p (O) is the closure of the smooth, compactly supported functions (j) : O —> R in 
the norm (J 0 \D(j)\ p dxj 1 ^ P ', we refer readers to the sources 01 for information on Sobolev 
spaces and their applications to PDE . It is evident that 1/A P is the smallest constant C for 
which the Poincare inequality 

I \i/j\ p dx <C I \Di(j\ p dx, tl> e 

Jn Jn 


holds. 
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The constant A p is also a type of eigenvalue. Indeed, minimizers in (II.ip are called 
p-ground states and satisfy the PDE 

{ —A p u = X p \u\ p ~ 2 u, x E Cl, 

u = 0, x E dVt. 

Here, the operator A p if> := dw(\Dfi\ p ~ 2 Dif) is known as the p-Laplacian. It has been 
established that p-ground states exist and that any two are multiples of one another, see 
mm- Consequently, A p is said to be simple. 

We will use the following iteration scheme to approximate X p and p-ground states. Let 
«o ^ L p (fl), and consider the family of PDE 

\—ApUk = \uk-i\ p ~ 2 Uk-i, xEfl . . 

Uk = 0, x E dQ 


for k E N . It can be verified without too much difficulty that for a given uq, there is a unique 
weak solution sequence (uk)keN C H / 0 1,p (O) of (jl.2j) . That is, there is a unique sequence 
(wfc)fceN C Wl' p (XV) such that 


\Duk\ p 2 Duk ■ D(j)dx = / \uk-i\ p 2 Uk-i4>dx 


for each 0 E W ( J’ P (Q) and k E N. In fact, once Uk-i E L P (Q) is known, Uk can be obtained 
by minimizing the functional 


1E 0 1,P (D) 3 v i y f f-\Dv\ p — |M fc -i| p V- iv) 


As this functional is strictly convex and coercive, the existence of a unique minimizer follows 
from the “direct method” of the calculus of variations. 

The following theorem details how the scheme (II.2p is related to A p and p-ground states. 

Theorem 1.1. Assume Uq E L p {kl ) and define 


Then the limit 



exists in Wq P (Q.). If ^ 0, then if is a p-ground state and 


Xp 


J n \Du k \ p dx 
k ^°° In \ u k\ p dx 


(1.3) 
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Remark 1.2. It may not be obvious how to verify that the limiting function ip is not identically 
zero. However, if for instance u 0 > 0 in H or if u 0 ^ 0 an d is regular enough in order to 
have a Hopf’s lemma (for instance C 1,a , cf. [ID]), then it is straightforward to verify that ip 
is indeed non-zero. 

The iteration scheme (11.2p was introduced by R. Biezuner, G. Ercole, and E. Martins in 
P] who conjectured the limit 


( J n \u k -i\ p dx \ 1 1/p 

P V In KN* J 


(1.4) 


We prove this limit holds under the hypotheses of Theorem 11.11 see Corollary 12.31 We also 
show that the sequences 


( and / f n 

V Jn\ u k\ p dx y fceN V In\ u k\ p dx / fceN 


are nonincreasing, which we regard as special features of the the iteration (11.21) . See Propo¬ 
sition 12.41 below. 

Next, we derive an iteration scheme in the limit as p —> oo. Our motivation was the 
seminal work of P. Juutinen, P. Lindqvist, and J. Manfredi |.6|, where it was proven that 
linip^oo A l J p exists and equals 


Aoo := inf 


\Dlp\L°°(fl) 

IV’I L°°(Cl) 


ipew!°°(n),ip^o 


(sup{r : B r (x) C for some x G O}) 1 


Here H / 0 1,oo (O) is the space of Lipschitz continuous functions ip : Q —* R that satisfy ip\on = 0. 
Furthermore, these authors also showed that there is a sequence (u Pj )j G ^ of p-ground states 
that converge uniformly to a viscosity solution w G H / 0 1,oo (O) of the PDE 


{ min{—AooW, \Dw\ — Aootc}, w > 0, 

-AooW, w = 0, 

max{—A oqW,—\D w\ — Xoow} : w<0. 


(1.5) 


Here ip := D 2 ipDip ■ Dip is the infinity Laplacian and nontrivial solutions of (ll.5[) having 
constant sign, are called oo-ground states. 

Passing to the limit as p —> oo in (11.21) . we are able to conclude the subsequent result. 
The novelty in the theorem below is that (j 1.6 [) presents a new mechanism for generating 
oo-ground states. 

Theorem 1.3. Assume Uq G (7(11) and denote (uk, p ) a,sn as the solution sequence of (II.2j) . 
(i) There is a sequence (pj)je^ increasing to oo and (vk)keN C W 0 1,oo (fi) such that —y Vk 

uniformly on fl as j —)■ oo for each fcGN. Moreover, Vk is a viscosity solution of the PDE 


{ min{—AooVfc, \Dv k \ - v k _ i > 0 

AooV k , v k —\ 0 

maxj-AooUfc, -\Dv k \ - v k -i}, ffc~i < 0 


( 1 . 6 ) 
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for each fceN. (Here u 0 := u 0 .) 

(ii) The limit L := lim^oo X^DvklL^in) exists. If L > 0, 

. \DVk\L°°(Q) 

Aqq lim j j . 

fc-)-oo \V k \ Loo(O) 

and any uniformly convergent subsequence of (X^v^k&N converges to a solution of (jl.5j) . 

Remark 1.4. Obviously, if Uq > 0 and L > 0, then any uniformly convergent subsequence of 
{Kx, v k)k<m converges to an oo-ground state. 

We would especially like to thank Richard Tapia. After learning about our previous work 
[5] which employed a doubly nonlinear flow to approximate X p and ground states, Professor 
Tapia suggested that it may be possible to use inverse iteration to obtain similar results. As 
noted above, the authors R. Biezuner, G. Ercole, and E. Martins were the first to make this 
observation in [lj. Nevertheless, we believe this paper adds significantly to [T] and makes 
clear the connection between inverse iteration and p-ground states. 


2 Convergence of the scheme 

Before proving Theorem 11.11 we will first make an observation which illuminates how /i p 
enters the statement of the theorem. In particular, we will argue that (n p Uk)keN is bounded 
in W 0 1,p (fi) and {^ p \Du k \LP(n)) k£N is a nonincreasing sequence of real numbers. 

Lemma 2.1. For each k E N, 

1 i p I \Du k+ i\ p dx < I \Du k \ p dx. 

Jn Jn 

Proof. Assume \ Du k +i\ p dx ^ 0. We employ Holder’s inequality and the Poincare inequal¬ 
ity to find 



Du k +i\ p dx 



\Du k+ i\ p 2 Duk+iDuk+idx 


< 


< 


\u k \ p 2 u k u k+1 dx 
_ i-i /p 

\u k \ p dx 


\ i /p 

\u k+ i\ p dx\ 


| Du k \ p dx 


•p Jn 


— ( / \Du k \ p dx 


i-i /p 


i-i /p 


- [ \Du k+1 \ p dx\ 
•p J u J 

\Du k+1 \ p dx\ 
n ) 


( 2 . 1 ) 
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Consequently, 


\Du k+ i\ p dx < 


\p/(p- 1) 
Ap 


\Duk\ p dx 


which proves the claim. 

Remark 2.2. A minor variation in the proof of Lemma [2.11 gives the estimate 

/ \Du k \ p dx < — f \uk-i\ p dx 
Jci Rp Jci 

for each k G N. This estimate will be employed in the proof of Theorem 11.31 
Proof of Theorem \l.l[ Set if k := R p Uk (k G N) and 


□ 


( 2 . 2 ) 


S := lim 

k—t OO 


\Dif k \ p dx. 


Observe that the limit defining S exists by Lemma 12.11 If S — 0, the assertion follows. Let 
us now assume otherwise. 

Notice that (if k )ken satisfies the sequence of PDE 

-ApVtfc = X p \if k -i\ p ~ 2 ifk-i, x G O, 

if k = 0, x g dd. 

By Lemma [All and Rellich-Kondrachov compactness, there is if G W / 0 ' :P (O) and a subsequence 
(ifkj so that if kj —>• if in L p (kl) and Dif kj Dif in L p (kl\ M”), as j — y oo. Also note 

I \Dif kj \ p dx = f \Dif kj \ p ~ 2 Dif kj ■ Difkjdx — X p f \'if kj - 1 \ p ~ 2 'ip kj -i'ip kj dx. 

Jci Jci Jci 

Since if kj —> if in L p (h2), 

lim sup f \Dif kj \ p dx = X p f \if\ p dx < f \Dif\ p dx. 
j->°o Jd 3 Jn 

And the weak convergence Dif kj Dif in L p (f2;M n ) gives 

| Dif k .\ p dx> f \Dif\ p dx. 

Jci 

Thus, ifkj —> if in S = f n | Dif\ p dx and 


lim inf 

Jci 


/ \Dif\ p dx = X p / | if\ p dx. 

J Cl J Cl 

As S > 0, if ^ 0 and thus if is a p-ground state. Since S is the same for all any subse- 
quential limit, the simplicity of \ p implies that if k if in W /7 Q 1 ’ P (0) as claimed. Moreover, 

J n \Du k \ p dx __ f Q \Dif k \ p dx _ f n \Dif\ p dx _ 

In \ u k\ p dx f n \ifk\ p dx f n \if\ p dx 

□ 
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Observe that if u 0 is a p-ground state, then (/i“ fc -u 0 )fcGN is a “separation of variables” 
solution of (11.21b This is a trivial case of Theorem 11.11 Also note that lim k^ooUpUk could 
vanish identically. For instance, this occurs when p = 2 and u$ is an eigenfunction of the 
Dirichlet Laplacian corresponding to an eigenvalue different that A 2 . Let us now see how the 
limit (II.4p follows from Theorem 11.11 

Corollary 2.3. Assume lim^oo Pp\Du k \LP(n) ^ 0, then the limit (II.4[) holds. 

Proof. Set ifk '■= hp u k- By the previous assertion, (VtOfceN converges to a p-ground state in 
VF 0 1,p (n). As a result, 


fnlu^dx = J n \'i/j k . 1 \Pdx 

°° LI u k \Pdx 'S-oo J n \A\ p dx 


\p/(p- 1) 
p 


□ 


We conclude this section by establishing some fundamental properties of the iteration 
scheme (11.21) . The monotonicity (12.31) suggests the iteration scheme improves the Rayleigh 
quotient J n \Dif\ p dx/ f Q [fj^dx at each step, and the monotonicity (12.41) gives more insight 
on the limit (II.4j) . 

Proposition 2.4. Assume uq G W 0 1,p (fi) and uq ^ 0. Then Uk ^ 0 for each k eN, 


L \Du k+1 \ p dx 
L \u k+1 \ p dx 

and 

L \u k \ p dx 
fa Wk+i\ p dx 

for each k E N. 


L \Du k \ p dx 

L \u k \ p dx 

(2.3) 

L \u k ~i\ p dx 

L 

(2.4) 


Proof. If Uo ^ 0, then u\ ^ 0 or (11.21) could not hold when k — 1. By induction, we may 
conclude u k ^ 0 for each k E N. 

Now £x k E N and observe 


\u k \ p dx = / (\u k \ p 2 u k )u k dx 


= / \Du k+l \ p 2 Du k+ i ■ Du k dx 
Jn 

/ r \ i- 1 /? / r 


i/p 


< 


\Du k+ i\ p dx 


| Du k \ p dx 


(2.5) 
























Combining the bound (12.ip with (12. 5 p gives 


fn \Du k +i\ p dx < (fgluklPdx ) 1 1/p (f Q \u k+ 1 \ p dx) L/p 
fn \uk+i\ p dx ~ J n \u k+1 \Pdx 

_ _ f Q \u k \ p dx _ 

(fn \u k +i\Pdx) l ~ 1/p (f n \u k \Pdx) l/p 
< (fn \Du k+ 1 \Pdx) 1 - l/p (f Q \Du k \Pdx) 1/p 

(fn K+iNx) 1-1/p (f n \u k \Pdx) 1/p 
= ( fn \Du k+ i\ p dx \ l ~ l/p / J^\Du k \Pdx \ 1,p 
V In \ u k+i\ p dx J V In \ u k\ p dx J 

which verifies (I2.3[) . 

As for (12.41) . we employ (12.51) . (12.3j) and (12.11) to find 


fn Wk\ p dx < 
fn \u k+ i\ p dx ~ 

< 


< 


(fn iDuk+^dx) 1 1,P (f Q \Du k \ p dx) 1/p 
Jn\uk+i\ p dx 


IJn 


\u k+1 \ p dx 

In 


f n \Du k \ p dx~ 
fn \ u k\ p dx . 
\Du k \ p dx 


i-i /p 


(In \Du k \ p dx) l/p 

fn \u k+ i\ p dx 


(In Wk+i\ p dxY /p (f n \u k \Pdx)^ 1/p 

(fn \u k \ p dx) 1/p (f Q \u k _i\ p dx) l l/p 

(In \uk+i\ p dx) 1/p (f Q \u k \Pdx) l ~ l/p 
f f n \u k \Pdx A 1/p / f Q \u k _i\ p dx \ 1 ~ 1/p 

\fn\ u k+i\ p dxj V f n \ u k\ p dx ) 


□ 


Remark 2.5. If u 0 ^ 0, the sequences 

( and / 

V fn \u k \ p dx J fcgN \ In \ u k\ p dx ) fcgN 

are bounded below by A p and respectively; see Proposition 2.8 of [Q. In view of the 

monotonicity (j2.3[) and (12.4p . both of these sequences are convergent. However, the limits 
(11.3 j) and (jl.4j) may not hold if lim*.-**, n p u k = 0. For example, these limits fail if p — 2 and 
Uq is an eigenfunction of the Dirichlet Laplacian that corresponds to an eigenvalue not equal 
to A 2 . 
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3 The large p limit 

This section is dedicated to a proof of Theorem 11.31 which characterizes the large p limit 
of the solutions of the iteration scheme (II.2ft . We begin with an important observation 
regarding weak solution sequences (upkeN C W 0 1,p (fi) of (11.21) when u 0 G C(fi) 

Lemma 3.1. Suppose u 0 G C(Q), and let (upken C Wl' p (Vl) denote the associated solution 
sequence of (11.2)1 . Then for each k G N, there is a.k G (0,1) such that 

«,G^(H)nr(n). 

Proof. It suffices to verify the claim for k = 1; the case k > 2 then follows from induction. 
Recall that (11.2)1 implies U\ G W 0 1,p (fi) is a weak solution of solution of 

\-A p ui = |mo| p_ 2 m 0 , iGfl, 

U\ = 0, x G dfl. 

We will use a weak comparison principle argument to bound u\ from above and then from 
below. The regularity theory developed by E. DiBenedetto in [3] would then imply the 
existence of an aq G (0,1) such that u\ G Cjjj” 1 (Q). 

To this end, we fix any y ^ 12 and define 

w(x) —\—\x — y\ g , x G fh 

qn~ 

Here q = p/{p — 1) is the Holder exponent dual to p. Direct computation has A p w(x) = 1 
for each iGO. It is also routine to verify that 

v ■= |wo|L°°(r2) (|'HL°°(n) — w ) 

satisfies 

— A p v > \u 0 \ p ~ 2 u 0 , x G H. 

Since u|an > 0 = ui|an, a standard weak comparison argument implies u\ < v in fl In 
particular, 

Ui < x G fl. 

We can argue similarly to bound u from below and derive 

Ml > — |rc|L°°(f 2 )|'Wo|L°°(r 2 ), x G fl 


□ 

We have just established that the solution sequence (upkm of the inverse iteration scheme 
is continuous, provided that u 0 is continuous. Virtually the same argument given by P. 
Juutinen, P. Lindqvist and J. Manfredi in the proof of Theorem 2.5 of [7] implies that each 












Uk is additionally a viscosity solution of (II.2p . That is, each solution sequence (w/JfceN C C(fi) 


of (ll.2p with p > 2 has the following property. For each fceN, 

-A p 0(x o ) < |'Ufc_i(x 0 )| P_2 Mfc-l(Xo) 

whenever <f G C' 2 (h2) and u k — f> has a local maximum at Xo G fl, and 

- A p (fi(x 0 ) > |rtfc_i(a;o)| p_2 Mfc-i(a;o) 

whenever f> G C' 2 (h2) and u k — <f> has a local minimum at £ 0 G hi. We refer interested readers 
to the “user’s guide to viscosity solutions” [2J for more information on viscosity solutions of 
elliptic PDE, and we are now ready to prove Theorem 11.31 

Proof of Theorem 1 1. 71 part ( i ). Employing Lemma [2. II and inequality ( 12.211 for k = 1 gives 



Assume po > n. For p > po, we can use the above inequality with Holder’s inequality to get 



By Morrey’s inequality and lim^oo/i p = Aoo, 

C C 1 _n/po (fi) 



is bounded for each fcGN. Therefore, the Arzela-Ascoli Theorem and a typical diagonaliza- 
tion argument implies there is a sequence (ufc)fceN C C 1 " n//po (hl) and a sequence of positive 
numbers ( pf)j ^n that are increasing and unbounded such that 


v k = lim u k . v . 


in C 1 n / po (hl) for each fcGN. 


Now let p > r, and employ Holder’s inequality and (12. 2 p to get 



- nwl' Ufc - 1 -pl L °°( f2 )- 
/V 
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The sequence ( u k, Pj )j>j r is then bounded in W 0 1,r (fi) for some j r G N large enough and thus 
converges to Vk weakly. Therefore, we can substitute p = p 3 above and send j oo to arrive 
at 

(pi J \Dv k \ r dx^j < |n fc _i| L oo(n). 

for each k G N. And after sending r —> oo, 


\Dvk\L°°(n) < |11• (3.1) 

In particular, we have verified that (vk)keN C W 0 1,oo (fi). 

We will now verify that Vk are viscosity solutions of the iteration scheme (11.61b By 
induction, it suffices to prove this for k = 1. Assume cp G C 2 (fi) and v\ — (p has a local 
maximum at Xq G fl. We aim to show, 


{ min{—A oo 0(x o ), \D(j)(x 0 )\ - Mo(^o)}, 

-Aoo0(xo), 

max{—A oo 0(x o ), -|D0(x o )| - n 0 (^o)}, 


uo(zo) > 0, 
uo(zo) = 0, 
u 0 (x 0 ) < 0. 


(3.2) 


After adding x > ||x — x 0 | 2 to 0 and later sending p —> 0 + , we may assume that v± — (p 
has a strict local maximum. Since U\ yPj converges to V\ uniformly on hi, there is a sequence 
(xj)jeN c converging to xq for which u\ iPj — (p has a local maximum at x 3 . Since ui tPj is a 
viscosity solution of ( 11 . 21 ) with k — 1 and p = Pj, 

-A Pj (p(xj) < \u 0 (x j )\ Pj ~ 2 u 0 (x j ). (3.3) 


If u 0 (x 0 ) < 0) then u 0 (xj ) < 0 for all j sufficiently large. By (13.3[) . 

- A Pj (p(x 3 ) = \Dcj)(x j )\ Pj ~ 4 {iD^Xj^A^Xj) + (pj - 2)A oa (p(x j )} < 0, (3.4) 

and thus \D(p{x 3 )\ Y 0 all large enough j G N. Canceling the factor of \D(p{x 3 )\ P:i ~ A in (13.4ft . 
dividing by pj — 2 and sending j —> 00 gives —A^cp^xo) < 0. Likewise, rearranging (13.31) 
leads to 


\D(p(xj)\ 2 A(j)(xj 


A 00 cp(x j ) < 


( \M x j) I 


Pj-A 


Uo(XjY 


(3.5) 


Pj ~ 2 — p j -2\\D(p{x j )\ 

Therefore, it must also be that —Uq(x 3 ) < \D(p(x 3 )\ for all j large enough. Hence, (j3.2[) holds 
when Uo(xo) < 0. 

Now suppose u 0 (x 0 ) = 0. If additionally, \Dcp(x 0 )\ = 0, then clearly — A oo 0(x o ) < 0. If 
|-D0(x o )| Y 0; we can send j —> 00 in (13. 5 p to again arrive at —A^cpf^xo) < 0. Thus, (j3.2j) 
holds when Uq(xq) = 0. 

Finally, let us assume that Wo(^o) > 0, and that \D(p(xo)\ — Uq(xq) > 0. Then \D(p(xj)\ — 
Uo(xj ) > 0 for all j 


G N sufficiently large. Passing to the limit in (13. 5 p again gives 
—Aoo0(xo) < 0. In conclusion, (13.2p holds in the case ^ 0 (^ 0 ) > 0, as well. Therefore, 
we have verified that V\ is a viscosity subsolution of (11.61) . An argument that shows v\ is 
additionally a viscosity supersolution of (II.6p can be made similarly, so we leave the details 
to the reader. □ 
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Proof of Theorem \1.3\ part (ii). In view of (13.ip . 


\Dv k \L°°{p.) < |ffc-i|z,°°(n) < t— \Dvk-i\L°°(n)- 

Aoo 


Therefore, the sequence (A^|.Dufc|i,°°(n))fc€N is nonincreasing, and the limit 


L := lim A^| Dv k \L°°(ii) 

k—> oo 


exists. The inequality (13.ip also implies 

|Wfc|L°°(n) < T — \DVk\L^{p) < T - 
/ '00 / '00 

Consequently, (A^ 0 |ufc|i / oo(Q)) fceN is nonincreasing and the limit 

M : = lim A^JufclL^n) 


exists, as well. 

Observe A^|Z>u fe | L cx, ( Q) < Aoo (A^I'Ufc-ili^n)) so that 

L < Aoo M. 

Moreover, A^l^l^n) < y^A^|T)u fe | L c» (n) , which implies 

Aoo M < L. 


Thus, Aoo M = L, and when this quantity is nonzero, 


Aoo 


lim 

k —^oo 


Dv k \L°°(Q) 

|Wfc|L°°(Q) 


Finally, note that that the sequence (w k ) ke ^ (A^oUfc)fcgN C W 0 1,oo (n) satisfies the 
iteration scheme 


0 


{ min{—AooWfc, \Dw k \ - AooWfc-i}, w k - 1 > 0 

Aqo w k , w k —\ 0 

max{—AooWfc, ~\Dw k \ - X^w^i}, w k _ x < 0 


in the sense of viscosity solutions. Therefore, if a subsequence of (X^Vkjk^ converges uni¬ 
formly on 0, the stability of viscosity solutions implies that the limit function is necessarily 
a solution of (II. 5p . □ 
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